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1.  INTRODUCTION 

The  phenomena  of  opto-acoustic  (OA)  pulse  generation  due  to  thermal  expansion 
induced  by  a  pulsed  optical  beam  has  been  known  for  some  twenty  years  now.  Earlier 
observations  of  OA  pulses  produced  by  pulsed  laser  beams  in  liquids  were  reported,  far 
example,  by  Askaryan  et  al.}  Carome  et  al ,,2  and  Rentzepis  and  Pao.3  Important  modern 
applications  of  pulsed  OA  generation  include  spectroscopy,4  energy  transfer  studies  and 
ultrasonic  measurements.5  However,  most  experimental  detection  of  the  OA  pulses  have 
been  performed  by  using  piezoelectric  transducers  or  microphones,  and  the  true  OA  pulse 
profile  produced  by  an  excitation  laser  pulse  cannot  be  readily  measured.  The  present  work 
seems  to  be  the  first  quantitative  experimental  investigation  of  the  true  OA  profile  (detected 
by  a  probe  beam  deflection  technique6,7,8)  as  a  function  of  the  excitation  laser  pulse 
parameters.  Such  quantitative  understanding  is  important  for  the  optimization  of  various  OA 
applications.5,9 


The  importance  of  theoretical  understanding  of  the  pulsed  thermo-elastic  OA 
generation  phenomenon  can  be  appreciated  from  the  many  theoretical  papers  on  this  topic, 
including  White,10  Goumay,11  Hu,12  Lyamshev  and  Naugol'nykh13  and  others.  Recently, 
explicit  mathematical  forms  for  the  OA  pulse  profile  in  a  liquid  excited  by  a  Gaussian  (spatial 
and  temporal)  laser  pulse  have  been  given  independently  by  Lai  and  Young14  and  Heritier.15 
Both  theoretical  treatments  have  provided  insight  as  to  the  method  to  produce  OA  pulses  of 
desired  profile;  for  example,  laser  pulses  of  narrow  diameter  and  short  duration  should 
produce  narrow  OA  pulses  propagating  in  the  lateral  direction  in  a  weakly  absorbing  medium. 
Their  theoretical  results,  which  turn  out  to  be  basically  equivalent,  are  described  in 
Section  2.  Our  experimental  arrangement  to  verify  these  theoretical  thermo-elastic  OA 
profiles  are  described  in  Section  3,  and  our  observations  are  presented  in  Section  4. 
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2.  THEORY 

Analytical  solution  of  the  pulsed  OA  generation  due  to  thermo-elastic  mechanism  in  a 
weakly  absorbing  medium  is  possible  for  the  simple  case  of  a  Gaussian  intensity  distribution 
in  space  and  time  of  the  excitation  laser  beam.  It  should  be  noted  that  the  exact  form  for 
the  Gaussian  distribution  may  vary  among  authors.  Lai  and  Young14  have  described  a 
Gaussian  laser  beam  in  the  following  form: 


(2*)3/2r\ 


exp 


'  2R2  2t?  / 


(1) 


where  ILY  is  the  laser  intensity  at  distance  r  from  the  beam  axis,  E  is  the  laser  pulse  energy, 
t  is  the  time  from  the  peak  of  the  laser  beam,  R  is  a  measure  of  the  laser  beam  radius  and  rp 
is  a  measure  of  the  pulse  duration.  Lai  and  Young  have  solved  the  acoustic  wave  equation 
with  a  source  term  being  due  to  the  thermal  expansion  produced  by  ILY(r.t),  and  obtained 
the  OA  pressure  Pi_Y(r,t)  in  terms  of  Bessel  function: 


PLY*r,t) 
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In  Eq.  (2),  K  determines  the  magnitude  of  the  OA  pulse,  re  is  a  time  scale  factor  which 
determines  the  width  and  the  magnitude  of  the  OA  pulse,  and  d$0/d£  determines  the  shape 
of  the  pulse.  The  symbols  used  have  the  following  meanings:  a  is  the  optical  absorption 
coefficient,  j5  is  the  volume  expansivity,  Cp  is  the  specific  heat  at  constant  pressure,  va  is  the 
acoustic  velocity,  ra=*=R/va  is  an  acoustic  transit  time  in  the  source,  t;  is  the  retarded  time 
with  t'«0  being  at  the  moment  r/va  after  the  laser  peak,  £  is  a  normalized  retarded  time, 
K1/4  and  Ij/4  are  Bessel  functions  of  imaginary  arguments  and  0  is  the  unit  step  function. 


Independent  of  Lai  and  Young,14  Heritier15  has  also  solved  the  thermo-elastic  wave 
equation  for  excitation  by  a  Gaussian  laser  beam,  with  an  intensity  distribution  IH(r,t) 


defined  by 


lH(r,t)  *  1/V1  exp(~  ‘  4) 

ir3/2w2r  '  w„  tz  ' 


Here  wQ  is  known  as  the  beam  waist  16  and  is  equal  to  2R,  and  r  is  the  e-1  half  width  of  the 
laser  pulse,  related  to  the  full  width  at  half  maximum  TpwHM  ty 


tFWHM  ”  *-665  T  • 


The  solution  for  the  OA  pressure  PH(r,t)  obtained  by  Heritier  due  to  the  thermal  expansion 
produced  by  IH(r,t)  in  a  weakly  absorbing  medium  is 

PH(r,t)  -  KV3/2F(x) 


where 


(2*)3/2Cfl 


or- 
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,  -  ,1/2 

e  -  (r2  +  w2/2v2)  (11) 

x  -  (t-  X.)/e  -  t'/e  (12) 

and 

F«-[r(?W-  hr*2) 

-  2xr(|-)1F1(i;|;x2)]  exp(-  x2)  .  (13) 

Here,  analogous  to  the  Lai  and  Young's  solution  in  Eqs.  (2)-(6),  Kf  is  an  amplitude  factor,  e 
is  a  time  scale  factor  and  F  is  the  shape  function  expressed  in  terms  of  confluent  (or 
degenerate)  hypergeometric  function  jFjfaibjz).  The  symbols  not  previously  introduced  are 
x  which  is  the  retarded  dimensionless  time  equal  to  the  retarded  time  t'  divided  by  e,  and  T 
which  is  the  Gamma  function. 

We  have  made  comparison  of  the  numerical  results  of  Lai  and  Young's  solution 
(Eqs.  (2-6))  and  those  of  Heritier's  solution  (Eqs.  (9-13)).  Both  solutions  state  that  the 
shape  of  the  OA  pressure  pulse  is  dependent  on  time  t  only  through  the  retarded 
dimensionless  time  £=tVre  in  Lai  and  Young's  result  and  x»t'A  in  Heritier's  result.  Due  to 
their  different  definitions  of  the  Gaussian  spatial  and  temporal  dependence  of  the  laser  pulse, 
it  should  be  noted  that  the  two  "time  scale  parameters"  e  and  re  differ  by  a  factor  of  /2: 

e  -  >/2tt .  (14) 

We  have  found  that  the  solution  for  ute  imt  : pendent  OA  pressure  profiles  PLY(r>0  *nd 
PH(r,t)  given  in  Eqs.  (2)  and  (9),  respectively,  are  essentially  equivalent;  this  is  shown  in 
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Fig.  1,  where  both  theoretical  profiles  are  plotted  with  respect  to  the  dimensionless  times  £ 
and  x. 


To  be  specific,  we  shall  use  Heritier's1^  notation  from  here  on.  e  is  called  the  "time 
scale  parameter"  of  the  OA  pulse,  because  e  truly  controls  the  temporal  width  of  the  OA 
pulse.  For  example,  Heritier's1^  line  shape  (Fig.  1)  shows  that  the  separation  between  the 
peak  and  the  trough  of  the  OA  pressure  pulse  is  Ax =1.66,  corresponding  to  a  time 
separation  being: 


Atpk.u  =  1.66s  . 


The  same  numerical  result  is  obtained  by  analyzing  Lai  and  Young's  results. 


(15) 


3.  EXPERIMENT 

The  present  experiment  is  a  first  attempt  to  verify  the  theoretical  OA  pulse  profile 
nredicted  by  Lai  and  Young14  and  Heritier.15  The  schematics  of  our  experimental 
arrangement  is  shown  in  Fig.  2.  Two  tjmes  of  excitation  lasers  with  different  pulse  durations 
are  used.  Most  data  were  obtained  with  a  Moiectron  N2  pumped  dye  laser  providing  laser 
pulses  of  TpwHM*®  nsec,  and  pulse  energy  E=0.16  mJ  at  10  Hz  repetition  rate.  The  dye 
laser  wavelength  is  580  nm  with  a  bandwidth  of  5  nm.  The  laser  beam  is  spatially  filtered 
and  an  iris  is  used  to  provide  a  4  mm  diameter  collimated  beam,  which  is  directed  to  a  90° 
prism  and  a  focussing  lens  mounted  on  a  translation  stage  (see  Fig.  2).  This  excitation  laser 
pulse  is  thus  focussed  into  the  sample  cell.  Some  of  our  data  have  also  been  obtained^  with 
an  excitation  laser  of  longer  duration  Msec)-  To  change  the  beam  waist  wQ  of 

the  excitation  laser,  we  use  focussing  lenses  of  different  focal  lengths  f.  For  the  N2  pumped 
dye  laser,  values  for  f  used  are  60  mm,  120  mm  and  254  mm;  the  corresponding  values  of 
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wQ  are  measured  by  using  a  calibrated  pinhole  to  be  20,  49  and  74  fim  respectively.  For  the 
flashlamp  pumped  dye  laser,  f  is  80  mm,  and  the  waist  w.  is  approximately  10  pm. 


The  probe  laser  beam  from  a  1  mW  HeNe  laser  is  beam  expanded  to  a  10  mm 
diameter  collimated  beam,  which  is  focussed  in  the  cell  by  a  lens  of  135  mm  focal  length, 
producing  a  probe  beam  waist  of  about  9  pm.  The  direction  of  the  probe  beam  is  adjusted 
to  be  parallel  to  that  of  the  excitation  beam.  This  is  done  by  first  moving  the  translation 
stage  controlling  the  excitation  laser  beam  so  that  the  two  beams  are  made  to  be  coincident. 
Subsequently,  the  two  beams  are  displaced  (by  moving  the  translation  stage)  to  a  separation 
of  typically  5  mm  for  the  present  measurements.  The  glass  sample  cell  is  of  rectangular 
shape  (5x3x8  cm)  and  contains  methanol  doped  with  a  dye  (eg,  Rhodamine  640)  at  about 
23  °C.  The  observed  OA  signal  shape  is  independent  of  the  dye  concentration  if  the 
transmission  of  each  laser  beam  is  not  less  than  ~30%. 


The  transmitted  probe  beam  is  directed  onto  a  6328a  interference  filter  (to  block  out 
scattered  light  from  the  dye  laser)  and  a  fast  photodiode-preamplifier  assembly  (Analog 
Module  model  LNVA-O).  The  small  active  area  of  the  photodiode  (0.8  mm2)  measures  the 
intensity  of  only  a  small  portion  of  the  probe  beam;  its  position  in  the  probe  beam  (of 
diameter  about  4  mm  at  the  photodiode)  can  be  adjusted  by  a  translation  stage.  When  the 
OA  pulse  produced  by  the  excitation  beam  crosses  the  probe  beam,  a  transient  angular 
deflection17  <p  of  the  probe  is  produced: 


*(r,t)  .  _L  «  <W>- 


where  n0  is  the  normal  refractive  index  of  the  liquid,  n(r,t)  is  the  change  in  the  refractive 
index  which  is  proportional  to  the  OA  pressure  P(r,t),  and  l  is  the  interaction  length  of  the 
OA  pulse  with  the  probe  laser.  The  small  transient  probe  deflection  <p  causes  the  probe 
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beam  to  move  across  the  detection  photodiode.  If  the  photodiode  is  suitably  positioned  at  a 
wing  of  the  probe  beam  cross  section,  the  deflection  0  produces  an  optimum  intensity 
change;  however,  if  it  is  positioned  at  the  center  of  the  probe  beam,  0  produces  much 
smaller  intensity  changes.  This  is  explained  in  Fig.  3,  where  the  observed  signals  are  also 
shown.  The  photodiode  and  preamp  assembly  has  a  bandwidth  of  100  Hz-100  MHz.  Its 
output  is  displayed  on  a  75  MHz  bandwidth  oscilloscope  (Tektronix  7854  with  7A24  plug 
in).  The  digitizing  and  signal  averaging  capability  of  the  oscilloscope  is  too  slow  to  be  useful 
for  recording  the  OA  signal;  hence,  any  signal  averaging  is  presently  performed  by 
photographing  the  scope  display  with  multiple  exposures  (eg.,  for  50  laser  shots). 

The  observed  probe-beam  deflection  signal  S(t)  from  the  photo-diode  is  given  by  the 
following  equation: 

S(t)  -  Gyrd)L0  ,  (18) 

where  G  is  a  constant  depending  on  the  photo-diode  sensitivity  and  electronic  gain  of  the 
/ 

detection  system,  Ip(rd)  is  the  lateral  spatial  derivative  of  the  probe-beam  intensity 
distribution  at  the  photo-diode  position  rd  (with  the  active  area  of  the  photo-diode  being 
sufficiently  small),  and  L  is  the  "lever  arm"  of  the  probe  beam  (i.e.,  distance  from  the 
interaction  region  in  the  cell  to  the  photo-diode).  Combining  Eq.  (17)  and  Eq.  (18),  we 
have 

S(t)  cc  dp(r,t)/  3t  ,  (19) 

which  means  that  our  experimental  probe-beam  deflection  signal  is  a  measure  of  the 
lime-derivaiivt  of  the  OA  pulse  at  the  probe-beam  position  in  the  liquid.  To  compare  with 
theory,  we  should  either  compare  the  observed  signal  with  the  timederivative  of  the 
theoretical  OA  profile,  or  equivalently,  compare  the  time-integral  of  the  observed  signal  with 


.i.  t 
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the  theoretical  OA  profile.  We  shall  choose  the  former  method,  since  we  believe  that  it  is 
better  to  make  theoretical  comparisons  with  the  raw  data  rather  than  with  "processed"  data 
wherever  possible. 

4.  EXPERIMENTAL  RESULTS 

The  experimentally  observed  probe  deflection  signals  for  three  different  focussing  of 
the  N2  pumped  dye  laser  are  shown  in  Fig.  4.  These  signals  are  obtained  with  tiie 
photodiode  being  in  a  wing  of  the  probe  beam  cross  section  (Fig.  3a)  so  that  the  signals  are 
maximized.  These  observed  signals  are  replotted  in  Fig.  5  in  dimensionless  times;  the  best 
fits  to  the  theoretical  deflection  signal  are  obtained  with  the  time  scale  parameter  e  being 
10.6,  30.1  and  42.9  nsec  for  focussing  lens  of  focal  length  60,  120  and  2S4  mm, 
respectively.  Our  results  are  summarized  in  Table  I,  which  gives  the  measured  beam  waist 
w0  and  the  corresponding  theoretical  values  of  e  according  to  Eq.  (11),  and  the  best  fit 
values  of  e.  We  see  that  the  theoretical  time  scale  parameter  according  to  Eq.  (11)  agrees 
with  the  value  obtained  by  fitting  the  experimental  profile  to  better  than  about  ±20%.  This 
provides  a  first  experimental  verification  of  the  theoretical  profiles  calculated.14*1* 

The  observed  probe  deflection  signal  for  the  case  of  flashlamp  pumped  dye  laser 
excitation  is  shown  in  Fig.  6;  the  best  fit  to  the  theoretical  signal  (Fig.  7)  is  obtained  with 
e*1.05  jusec.  Here,  if  we  fit  the  excitation  laser  pulse  of  tfwhM*1-3  (see  Fig.  6)  by 
a  Gaussian  intensity  distribution,  we  obtain  a  laser  width  parameter  t«780  nsec.  For  such  a 
long  pulse  duration  and  small  beam  waist  (wo«10  /am),  Eq.  (11)  indicates  that  e  should  be 
approximately  equal  to  r.  The  ~2 5%  difference  between  the  fitted  OA  time  scale  parameter 
e  and  laser  pulse  width  r  is  partially  because  the  laser  intensity  is  significantly  nonGaussian. 
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It  is  clear  from  Figs.  S  and  7  that  there  is  an  imperfect  agreement  between  the 
experimentally  observed  OA  signal  and  the  theoretical  signal  predicted  by  Lai  and  Young14 
or  by  Heritier.15  The  agreement  is  better  for  the  earlier  half  of  the  signal  than  the  latter 
half.  The  time  scale  parameter  e  obtained  by  fitting  the  observed  signal  to  the  theoretical 
one  agrees  with  the  value  calculated  from  Eq.  (11)  to  better  than  ±30%.  Some  of  the 
reasons  for  the  limited  agreement  are  the  following:  (a)  Neither  the  spatial  nor  the  temporal 
dependence  of  the  excitation  laser  pulse  is  exactly  Gaussian,  as  the  theoretical  results 
assumed;  (b)  The  excitation  and  probe  laser  beams  may  not  be  exactly  parallel,  and  they  are 
also  not  of  infinite  lengths  with  uniform  cross  sections,  so  that  the  theoretical  one-dimension 
treatment  is  not  totally  valid. 

5.  CONCLUSION 

We  have  used  a  probe  beam  deflection  technique  to  measure  the  time-derivative  of  the 
pulsed  OA  signal  in  a  liquid  generated  by  an  excitation  laser  pulse  of  tfwhm  being  8  nsec 
or  1.3  fisec.  The  observed  signals  are  in  reasonable  agreement  with  the  theoretical  results  of 
Lai  and  Young14  and  Heritier,15  which  we  have  shown  numerically  to  be  equivalent.18  We 
have  thus  experimentally  confirmed  (to  an  accuracy  of  ±30%)  that  the  width  of  the  pulsed 
OA  signal  is  controlled  by  the  magnitude  of  re,  as  given  by  Eq.  (4),  or  equivalently  by  the 
magnitude  of  e,  as  given  by  Eq.  (11).  Thus,  sharp  OA  pulse  is  generated  by  a  narrow  laser 
pulse  of  short  duration  r  and  focussed  to  a  small  beam  waist  wQ.  For  example,  conunerical 
lasers  are  readily  available  with  t«0.5  nsec;  if  the  laser  beam  is  focussed  to  a  waist  wQ  of 
1  fan,  this  would  produce  an  OA  signal  with  a  time  scale  parameter  e«0.8  nsec.  Since  this 
corresponds  to  a  peak  to  trough  time  separation  of  1.33  nsec  (see  Eq.  IS),  a  very  sharp 
acoustic  pulse  with  Fourier  frequencies  into  the  GHz  regime  is  produced  in  the  sample.  The 
monitoring  of  the  propagation  of  such  an  acoustic  pulse  In  the  sample  should  provide  a  novel 


method  of  frequency  multiplexed  measurement  of  ultrasonic  dispersion  and  attenuations  in 
matter.9 
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TABLE  I 

Summary  of  Experimental  Parameters  for  Our  Study  of  OA  Pulse  Generation 
by  a  N2  Pumped  Dye  Laser  Beam  (r— 4.8  nsec). 


Focal  length  f 

Measured  beam 
waist  wQ 

e  according  to 

Eq  (11)  with 
measured  wQ  and 
t*4.8  nsec 

e  obtained 
by  a  best 

theoretical  fit  of  the 
observed  probe 
deflection  signal 

60  (mm) 

20  (fim) 

13.6  (nsec) 

10.6  (nsec) 

120 

49 

32.3 

30.1 

254 

74 

47.6 

42.9 
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Figure  3.  Observed  probe  beam  deflection  signal  depends  on  the  position  of  the  small 
photodiode  in  the  probe  beam  cross  section.  This  is  shown  in  the  observed  signals  shown  in 
the  scope  photos  for  the  photodiode  being  on  one  wing,  on  center  and  on  the  other  wing  of 
the  probe  beam  cross  section,  respectively.  The  inset  explains  how  a  deflection  <j>  of  the 
probe  beam  produces  different  photodiode  signals  for  the  different  cases. 


Figure  5.  Comparison  of  theoretical  probe  beam  deflection  signal  with  the  observed  signals 
of  Fig.  4  plotted  against  the  dimensionless  time  x,  with  t  used  being  10.6,  30. 1  and 
42.9  nsec  for  (a),  (b)  and  (c),  respectively. 
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Figure  6.  Observed  probe  beam  deflection  signal  for  excitation  by  the  flashlamp  pumped  dye 
laser.  The  inset  shows  the  temporal  profile  of  the  excitation  laser  pulse. 


